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Abstract 

We investigate the thermal and magnetic properties of the integrable su{A) ladder model by means 
of the quantum transfer matrix method. The magnetic susceptibility, specific heat, magnetic 
entropy and high field magnetization are evaluated from the free energy derived via the recently 
proposed method of high temperature expansion for exactly solved models. We show that the 
integrable model can be used to describe the physics of the strong coupling ladder compounds. 
Excellent agreement is seen between the theoretical results and the experimental data for the 
known ladder compounds (5IAP)2CuBr4-2H20, Cu2(C5Hi2N2)2Cl4 etc. 

PACS numbers: 75.10.Jm,64.40.Cn 
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The experimental realization of compounds with a ladder-like structure has con- 
tributed to the intense interest in low-dimensional quantum systems. The existence of 
a spin gap, magnetization plateaux, quantum critical points and superconductivity under 
hole doping are examples of key physical properties observed in the ladder compounds. Of 
particular importance are the properties under a magnetic field H . According to the per- 
turbation theory result , the first-order terms for the zero temperature energy gap A 
and the critical field ifc2 are given in terms of the rung {J a.) and leg (Jy) exchange cou- 
plings by A = — J|| and fiBgHc2 = J± + 2J||. These results are in good agreement with 
the experimental data for strong coupling ladder compounds. However, the calculation of 
properties such as the full temperature phase diagram, the high field magnetization curve 
and the specific heat provide a significant challenge. 

We demonstrate here that the integrable su{4) ladder model 0,0] is capable of describing 
the physics of the ladder compounds. Indeed, the thermodynamic Bethe Ansatz (TBA) 
applied to the integrable sm(4) ladder model predicts the critical fields Hd = J± — 4J|j/7 and 
liBgHc2 = J±+4J||/7^here 7 is a rescaling parameter, which are also good fits for the strong 
coupling compounds p . Very recently the high temperature expansion (HTE) method 0, [J 
has suggested a way to calculate the full thermodynamic properties of integrable models 
Tom the so called T-system j8|] appearing in the quantum transfer matrix (QTM) formalism 
^, . Here we extend this approach to the integrable ladder model to derive the thermal 
and magnetic properties of the strong coupling ladder compounds. We compare our results 
with the experimental data obtained for the compounds Cu2(C5Hi2N2)2Cl4 
(5IAP)2CuBr4-2H20 H- 

The Hamiltonian of the ladder model is 
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H = — Hicg+J± ^ Sjfj—fiBgH'^^i^j+'^j)^ ^leg = ^ + + 4 SjSj+if, 

^ i=i i=i i=i 

(1) 

where Sj and Tj are Heisenberg operators, /i^ is the Bohr magneton and g is Lande factor. 
Throughout, L is the number of rungs and periodic boundary conditions are imposed. In the 
strong coupling limit, the contribution to the low temperature physics from the multi-body 
term in Tiieg is minimal and, as a consequence, the integrable ladder Hamiltonian exhibits 
similar critical behavior to the standard Heisenberg ladder . We adapt the model into the 
QTM method ^. The eigenvalue of the QTM (up to a constant) is obtained by the nested 



Bethe Ansatz to be 
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In this equation the chemical potential terms are /ii = Jx/2, ^2 = I^bqH, /X3 = and 
/i4 = —fiBgH, with = {v ± iu^)^ ■ The inhomogeniety parameter mat = — with 



'^al'^) = ni!li '('^"li"'')) = 1,2,3. Here N denotes the Trotter-Suzuki number. The 



fused Tm^ system j8|, which denotes the row-to-row transfer matrix with fusion type (a,m) 
in the auxiliary space carrying the m-fold symmetric tensor of the a-th fundamental repre- 
sentation of the 5^(4) algebra is essentially generated by the QTM eigenvalue T^^^ in Q. 
Thus rf'^^ can be embedded into the fused Tm ■* system. The analytic non-zero and constant 
asymptotic properties of the normalized Tm\v) system suggest the expansion ansatz 
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^^^iv) = ^^"o^ c^nW^/(i^-^ + {a + 1)2/4)". The QTM eigenvalue satisfies a set of the 



nonlinear integral equations 
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Following Ref. [3], the coefficients c)^', can be obtained recursively from Eq. with initial 



conditions 60"^ = InQ^"'', where Q^^' are constants related to the chemical potential terms 



via lim 



gS"^ with Q^^' = 1 and Qf' = exp{J^/2T). In this way the 
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spin ladder free energy /(T, //) 
first few orders we have 



can be expanded in powers of J\\/^T. For the 
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where A = 3,^(1 + 2Bo,i)/B'^^ with e = exp(Jx/2T) and 
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We find that the analytic expression (0) is sufficiently accurate to evaluate the model's 
thermodynamics. Nevertheless, we have considered the HTE up to fifth order. 

The experimental measurements of the susceptibility and magnetization of the compound 
(5IAP)2CuBr4-2H20 (abbreviated B5i2aT) suggest a spin ladder with exchange cou- 
plings = 13. OK and Jy = 1.15K. From the HTE for the integrable model we find that 
the values J± = 13. 3K and J\\ = 1.15K with rescaling parameter 7 = 4 give excellent fits 
to both the susceptibility and magnetization The temperature dependence of the sus- 
ceptibility is shown in Fig. ^ The agreement with the theoretical curve derived from the 
HTE is clearly excellent. The typical rounded peak for low magnetic field, characteristic of 
a low dimensional antiferromagnet, is observed around 8. IK, in excellent agreement with 
the experimentally estimated value of 8K. The inset in Fig. ^ shows the low temperature 
behavior of the susceptibility, which is in excellent agreement with the experimental data. 



The theoretical curves for the high field magnetization shown in Fig. |21for different tem- 
peratures are also in good agreement with the experimental values. The field dependent 
magnetization curve predicts the low temperature phase diagram as well as the magneti- 
zation plateaux. For very low temperature the rung singlet forms a dimerized groundstate 
if the magnetic field is below the critical field Hd- The length of the antiferromagnetic 
correlation is finite while the triplet state is gapfuU. For finite temperatures the triplet 
excitations are also involved in the gapped phase. This can be observed in the high field 
magnetization curves for T = 1.59K and T = 4.35K in Fig. |21 At the critical field Hd, the 
gap is closed with fisgHd = A. If the magnetic field is above the critical point Hd, the 
lower triplet component becomes involved in the groundstate. For zero temperature, it can 
be rigorously shown that the other two higher triplet components do not become involved 
in the groundstate i5|. It follows that the strong coupling ladder can be mapped to the 
XX Z Heisenberg chain with an effective magnetic field term ^ Q|. The magnetization 
increases almost linearly with the field towards the critical point Hc2, where the ground 
state is fully polarized. At T = 0.4K, the HTE magnetization curve indicates Hd ~ 8.3T 
and Hc2 ~ 10. 5T, which are in excellent agreement with the experimental estimates of 8.4T 
and 10. 4T. The experimental magnetization in the singlet groundstate at low temperature 
appears to be nonzero. This nonzero magnetization attributed to paramagnetic impurities, 
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FIG. 1: The susceptibility versus temperature for B5i2aT at i/ = IT UM. The solid line denotes 
the susceptibility evaluated directly from the HTE. A parameter fit suggests the coupling constants 
J_L = 13. 3K and Jy = 1.15K with 7 = 4, = 2.1 and hb = 0.672K/T. The inset shows the same 
fit to the susceptibility at low temperature. 



which drive the low-temperature deviation between the experimental and the theoretical 
curves in Fig. El 

The inflection point is clearly visible in the experimental magnetization curves . This 
point is also evident in the theoretical curves at fisgH ~ J± where the magnetization 
moment is ^. The physical meaning of the inflection point is that the probabilities of 
the singlet and the triplet states | t|) in the groundstate are equal. Therefore, for the 
strong coupling ladder compounds at zero temperature the one-point-correlation function 



lies in a gapped singlet groundstate, which indicates an ordered dimer phase. 



while {Sj ■ Tj) = | in the fully-polarized ferromagnetic phase. However, in a Luttinger 
liquid phase, we flnd {Sj ■ Tj) = — | + S^. At low temperatures T -C J±, the one point 



correlation function is given by {Sj ■ Tj) 
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phase transitions can be clearly seen from the one-point correlation function curve shown in 
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FIG. 2: Magnetization versus magnetic field for B5i2aT with the same constants as in Fig. ^ 
The discrepancy in the magnetization curves at T = 0.4K and T = 1.59K is due to paramag- 
netic impurities which become negligible for higher temperature. The inset shows the one-point 
correlation function vs magnetic field. 



the inset of Fig. El 



We turn now to the ladder compound Cu2(C5Hi2N2)2 CI4 |ll|, |l2| (abbreviated 
Cu(Hp)Cl). In Fig. El we show the zero field magnetic susceptibility curve obtained from 
the fifth order HTE free energy. A full fit with the experimental data suggests the coupling 
constants J_l = 13. 5K, Jy = 2.4K and 7 = 5. The effect of the magnetic field is to lift the 
susceptibility in the low temperature regime. We notice that there is a discrepancy with the 
zero temperature TEA result for lic2 H due to the presence of strong exchange coupling 
along the legs. For finite temperature this discrepancy is smaller. The inset curves in Fig. El 
show the high field magnetization for temperatures T = 0.6, 1.6 and 4.04K. We observe that 
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FIG. 3: Susceptibility versus temperature for the compound Cu(Hp)Cl Q]. The solid line denotes 
the susceptibility evaluated directly from the HTE with = 0.672K/T, J_l = 13. 5K, Jy = 
2.4K, 7 = 5 and g = 2.03. The inset shows the magnetization versus magnetic field at different 
temperature. At T = 0.6K, the critical fields are Hd ~ 7.8T and Hc2 ~ 13. OT, in good agreement 



with the experimental results 
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experimental values 
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7.8T and i7c2 ~ 13. OT, which are in good agreement with the 



12| . It is also obvious that the finite temperature causes a spin-flip 
in the gapped ground state. 

The specific heat curves in Fig. 0] for i7 = OT and H = 4T indicate that the HTE result 
also agrees satisfactorily with the experimental data 3| • In the absence of a magnetic field, 
a rounded peak indicating short range ordering is observed around 4.5K. For temperatures 
T < 4.5K, there is an exponential decay due to an ordered phase. The humps become 
smaller as the magnetic field increases. For the H = 4T curve a peak is observed at around 
4K. As to be expected, there appears to be a small deviation from the experimental data 
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FIG. 4: Specific Heat versus temperature for the compound Cu(Hp)Cl 12] with the same constants 
as in Fig. El The inset shows the field dependent entropy versus temperature. 

at very low temperatures. The inset of Fig. |3] shows that the entropy curves for magnetic 
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fields H = OT and H = 4T are also in agreement with the experimental data 

The full phase diagram for the two compounds is shown in Fig. El The slopes of the critical 
curves indicate that the estimated values of Hd and Hc2 at T = OK coincide with the TBA 
results j^. We have also examined other strong coupling compounds. Comparison with the 
experimental data for the compound BIP-BNO [l2| suggests the coupling constants J± = 
75K, J|| = 15K with 7 = 6.0 and ^ = 2.0. For the compound [Cu2(C202)(CioH8N2)2)](N03)2 
jisl ] we find J± = 515K, Jy = 40K with 7 = 5.0 and g = 2.14. The respective spin excitation 
gaps, A ^ 52K and A ^ 460K, are also in good agreement with the experimental values. 
The application of the HTE method to other ladder models, such as the mixed-spin ladders 
should now be straightforward. We also note that the HTE method predicts a fractional 
magnetization plateau with respect to different Lande (7-factors in the su{A) spin-orbital 
model, which coincides with the Bethe ansatz |3] and TBA 2(i| results. 
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FIG. 5: Phase diagrams for the compounds B5i2aT and Cu(Hp)Cl. 
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